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ABSTRACT 



We report new statistical equilibrium calculations for Fe I and Fe II in the 
atmosphere of Late- Type stars. We used atomic models for Fe I and Fe II having 
respectively 256 and 190 levels, as well as 2117 and 3443 radiative transitions. 
Photoionization cross-sections are from the Iron Project. These atomic models 
were used to investigate NLTE (non local themodynamic equilibrium) effects 
in iron abundances of Late- Type stars with different atmospheric parameters. 
We found that most Fe I lines in metal-poor stars are formed in conditions 
far from LTE (local thermodynamic equilibrium). We derived metallicity 
corrections of about 0.3 dex with respect to LTE values, for the case of stars 
with [Fe/H] ~ —3.0. Fe II is found not to be affected by significant NLTE 
effects. The main NLTE effect invoked in the case of Fe I is overionization by 
ultraviolet radiation, thus classical ionization equilibrium is far to be satisfied. 
An important consequence is that surface gravities derived by LTE analysis 
are in error and should be corrected before final abundances corrections. This 
apparently solves the observed discrepancy between spectroscopic surface 
gravities derived by LTE analyses and those derived from Hipparcos parallaxes. 
A table of NLTE [Fe/H] and logg values for a sample of metal-poor late-type 
stars is given. 

Subject headings: NLTE effects — abundances — surface gravities 
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1. Introduction 

Iron is a basic stone for the study of the chemical evolution of stellar systems. Relations 
between elemental abundance ratios [X/Fe] versus [Fe/H] Q of giant and dwarf stars are 
generally used as tracers of chemical evolution of galaxies. The reason of this choice is 
that Fe lines are often quite numerous and easy to detect, even in very metal-poor stars. 
Hence, a good and precise determination of Fe abundances is of fundamental importance. 
Our understanding of Fe abundance in stars is mainly based on local thermodynamic 
equilibrium (LTE) analyses, for which many weak Fe I and Fe II lines are used. Most of the 
works devoted to spectral abundance analysis assume that the majority of these weak lines 
are not affected by important LTE deviations, in particular in the case of solar type dwarf 
stars. However, empirical studies of the ionization equilibrium of Fe I show some evidences 
that important NLTE (non local themodynamic equilibrium) effects are suspected to be 
present in the photosphere of very metal-poor stars (Zhao & Magain 1990, Fuhrmann 1998, 
Feltzing & Gustafsson 1998). Calculations performed by Takeda (1991), which refers largely 
to other excellent computational works in the literature, imply in marginal NLTE effects 
on the ionization equilibrium for the giant star Arcturus and more important ones when 
removing the UV line opacities. NLTE effects were also investigated by Gigas (1986) in the 
hot star Vega, finding a correction of 0.3 dex for the iron abundance with respect to LTE 
value. Rutten (1988) and Steenbock & Holweger (1984) have also enlighten the problem, 
determining the most important NLTE mechanisms working in stellar atmospheres of 
Late- Type stars. If important corrections should have to be applied to LTE abundances, 
this could have a real importance on galactic chemical evolution models. 

The basic theoretical problem of determination of stellar abundances through high 
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resolution spectral analysis is to solve the equation of radiation transfer, giving the variation 
of the radiative energy flow throughout an absorbing and emitting gaseous medium. 
Firstly, the solution requires the construction of a stellar atmosphere model, which gives 
the thermodynamic variables, temperature T and pressure P, as a function of the optical 
depth r. Secondly, we have to determine the absorption and emission coefficients i^ and k^,, 
which are directly proportional to the transition probabilities and the number of atoms, 
ions or molecules in a given quantum state or energy level. Atomic, ionic and molecular 
populations depend on the elemental abundance and on the degree of excitation, ionization 
and dissociation respectively. These last quantities are calculated through the solution 
of the statistical equilibrium equations for given conditions of abundances, gas density, 
temperature, depending also on atomic and molecular constants. In a NLTE problem, 
elemental populations will depend, not only on gas temperature and density (as in LTE 
case), but also on the radiation field. To take this into account, we must analyze each line 
transition by considering the transitions of neighbouring levels and not as a two-level atom 
as in the LTE case. As a consequence, we have to calculate simultaneously the transfer and 
the statistical equihbrium equations for each considered level. 

The real difficulty to treat a complete iron atom in the NLTE case, is to solve the 
statistical equilibrium equations including about 300 terms and 5000 multiplets. The 
diagnostic concerns a lot of strong UV lines, a lot of optical lines (as is used in classical 
detailed analysis), and also a lot of infrared lines which come from the highest terms of 
the atomic model, i.e., those with an excitation potential greater than 3 eV. We present 
in this work a much more complex model for Fe I and Fe II, taking into account a larger 
number of levels and transitions than those found previously in the literature. Our main 
goal is to estimate NLTE effects in the determination of Fe abundances using the curve of 
growth technique, as in a classical LTE analysis. In addition, this procedure allow us to 
check the ionization equilibrium and, consequently, to evaluate NLTE effects in the logg 
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determination. 

In section 2 we present our atomic models and the strategy to correct Fe abundances, 
and in addition surface gravities, of NLTE effects. In sections 3 and 4 we discuss the results 
for the Sun and stars previously studied in other works, respectively. In section 5, our 
results for metal-poor stars are given, and in section 6 we present our conclusions and some 
suggestions for future stellar abundance analyses. 

2. The iron atom model and the strategy 

2.1. Fe I and Fe II atoms 

Our first step for the analysis of NLTE effects was the elaboration of Fe I and Fe II 
atomic models. Our goal was to construct the best models from the statistical point of 
view, taking into account a complete set of levels and transitions. We have not tried to 
produce synthetic spectra to compare to observed profiles, but to give differential abundance 
corrections using equivalent widths and curves of growth, as will be described later in 
section 3. 

For Fe I and Fe II we consider levels with principal quantum numbers n=l,3,5, 7 and 
n=2,4,6,8 respectively. For computing time reasons we restricted our Fe I /Fe II models to 
256/190 levels, having a potential less than 6.8/8.7 eV and a continuum at 7.90/16.1 eV. In 
both models, fine structure were taken into account, resulting in a total of 2117 radiative 
transitions for Fe I and 3443 for Fe II. We have not considered line transitions in the far 
infrared (A > 5.0fim). Grotrian diagrams of our Fe I and Fe II models are given in Fig. 1. 

The code used to solve the equations of statistical equilibrium and radiative transfer 
is MULTI, described by Carlson (1986), which use the operator perturbation technique 
of Scharmer & Carlsson (1985). We used version 2.2 (1995). Radiation fields for each 
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stellar atmospheric model are computed using opacities from Uppsala package, including 
effects of some line blanketing. We did not use the option of the code which produces 
NLTE background opacities, once the models of atmospheres used here are based on LTE 
calculations (see 2.2) and our results do not have such degree of precision for individual 
line profiles yet (see section 3). As stated above, we are interested in the derivation of 
differential abundance corrections based on a curves of growth analysis. A better treatment 
of the opacity would be its derivation with the same code used to compute the atmospheric 
model, however, such treatment would probably not change strongly our results, mainly for 
metal-poor stars. 

Input atomic data for this code are: 

a) Energy levels: 

• excitation potentials 

• statistical weights 

• ionization stages 

b) Transitions: 

• oscillator strengths 

• radiative and collisional damping coefficients 

• photoionization cross-sections 

• excitation and ionization collisional cross-sections. 

Excitation potentials of the levels and their statistical weights are given by Hirata 
& Horaguchi's table (1995). For radiative transitions we used oscillator strengths given 
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by Fuhr, Martin & Wiese (1988), Hirata & Horaguchi (1995) and Thevenin (1989, 1990). 
Damping coefficients are used for calculation of line profiles (line broadening). The total 
line damping coefficient is given by 7 = 7rad + 7coii, where 7rad and 7coii are respectively 
the radiative and collisional damping coefficients. Radiative or natural damping width 
is defined as : 7rad = Z]i<i(Aii) + Z]i<j(Aji) where Au are the Einstein coefficients. The 
collisional damping coefficient is the sum of van der Waals and Stark coefficients, which 
take into account effects due to perturbations with neutral H and He (van der Waals) and 
charged particles (Stark). For all lines the classical van der Waals damping is evaluated 
from the classical approximation (Unsold 1955) and multiplied by an enhancement factor, 
since van der Waals constant cannot reproduce the real lines profiles (see Gurtovenko & 
Kostic 1982, Thevenin 1989 or more recently Anstee, O'Mara & Ross 1997), which will be 
discussed in section 3. In our case, damping due to Stark effect can be neglected, because 
the electronic density in the photosphere of the Late- Type stars investigated here is much 
smaller than neutral H density. 

Photoionization for all levels of Fe I and Fe II was treated in details by using the 
frequency dependence of the cross-sections given by the Iron Project (Bautista 1997). At 
each layer of the star's photosphere models ionizing radiation fields are computed, giving an 
estimate of photoionization rates. Note that if the ionization equilibrium is not achieved, 
the opacity computed would have to be changed by changing the iron abundance. The very 
detailed cross-sections given by the Iron Project were smoothed to decrease the number of 
points in tables cross-section vs. frequency for each level, in order to reduce the computing 
time. It was easy to check on some levels that this smoothing of the detailed cross-sections 
had no consequences on the final results of the radiative transfer computations. But, for 
some of the levels, when strong resonances in the photoionization cross-sections near the 
threshold exist, details were considered. 



CoUisional ionization rates with electrons are derived from the approximate formulae 
given by Mihalas (1978). For the evaluation of collisional rates of permitted transition lines 
we used van Regemorter's formulae (1962). For some forbidden lines considered, collisional 
rates were derived from the formulae of Auer & Mihalas (1973) with ^2=0.1, not equal to 1.0 
as proposed by Takeda (3.1.2, 1991)(see Cayrel et al. 1996). However, these forbidden lines 
seem to play a negligible role in our Fe models. Uncertainties in the collisional processes 
with neutral H and He are probably the main source of errors in our models. The lack 
of accurate cross-sections for collisions with hydrogen atoms is largely discussed in the 
recent literature. The use of the Drawin's theory (1968, 1969), as proposed by Steenbock & 
Holweger (1984), was severely criticized by Severino, Caccin & Gomez (1993), estimating 
that Drawin's theory gives cross-section values larger by a factor of lO'^ for NaD lines. A 
clear sum-up of the situation is given by Holweger (1996). In the absence of a more precise 
theory we used this approach as already done in many other works. But, we emphasize 
that we found globally no important consequences on the population distribution in atomic 
levels and therefore on the equivalent width of lines transitions of our Fe 1 and Fe II models, 
when using or not hydrogen collisions. 

We notice that changing Bautista's photoionization cross-section values by a factor 2 
induces changes in the resulting populations, giving errors on [Fel/Fell] less than 0.02 dex. 
This means that our strategy developed around the technique of the curve of growth (see 
2.3) has a low dependence on such atomic parameters. Of course this is not true if profiles 
of Fe lines are used to perform the stellar chemical abundance analyses. 

2.2. The atmospheric models 

In order to investigate NLTE effects on the ionization balance of iron with different 
stellar atmospheric parameters, atmosphere models for different known stars were generated 
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using the Bell et al. (1976) grid. This grid was used to be self-consistent with Thevenin's 
catalogue (1998), which reanalyzes LTE abundances of 35 chemical elements for 1107 stars. 
These models reproduce well the atmosphere of Pop 11 F to K dwarfs and giants, for 
which LTE atmospheric parameters were taken from Thevenin's catalogue. Microturbulent 
velocities were taken independent of the optical depth. No macroturbulent velocity was 
used in our computations. 

We have check that Bell et al.'s models (1976) and Kurucz's models (1993) for a dwarf 
metal-poor star give no significant differences on results of [Fel/Fell]. 

The solar iron abundance adopted in this work is Ape = 7.46 (Holweger 1979) as also 
used by Thevenin (1998). The reason of this choice is clearly explained in Thevenin (1989, 
see discussion & Table 111). Other determinations have lead to similar solar abundances 
(see for example Holweger, Heise & Kock 1990 or Holweger et al. 1991). 

2.3. The strategy to estimate NLTE effects 

To compute NLTE effects on the derived LTE surface gravities and [Fe/H] values, we 
constructed NLTE curves of growth for both Fe 1 and Fe II, using the equivalent widths 
W calculated with MULTl for a given stellar atmospheric model and for all transitions 
considered, as a simulation of observed equivalent widths . The abscissae of the curves of 
growth are computed under the condition of LTE as classically done by detailed analysis, 
for each of the 2117 Fe 1 and 3445 Fe 11 lines. However, to fit the classical LTE curves of 
growth for each analyzed star, we used only lines ranging between A A 2200-10000 A to be 
consistent with classical detailed analysis procedure. Once that atmospheric parameters 
estimated for a given star are from LTE analysis (see Table 2, Thevenin 1998), we can 
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immediately check the vahdity of this assumption and correct, if necessary, LTE values of 
[Fe/H] in the same table. But, before to correct the metallicities, we have to check the 
ionization equilibrium to derive the error on log g due to NLTE effects. If [Fel/H] were not 
equal to [Fell/H], this would mean that the classical LTE detailed analysis is wrong and 
would have forced the ionization equilibrium using a wrong surface gravity. Hence, for all 
stars, we had to calculate theoretical curves of growth with new corrected values of logg 
and then estimate metallicity corrections. 

3. Results in the solar case 

For the Sun, we used a model from the grid of Bell et al. (1976) given by Gustafsson 
(1981) and, in addition, tested the Holweger-Miiller's model (1974) to compare the results 
of both models for profiles of strong lines. We used the same solar abundance 7.46 for 
both models. Defining Wnlte and Wlte as the computed equivalent width in NLTE and 
LTE conditions respectively, results of Wnlte/Wlte ratios, which give the importance of 
NLTE effects on each spectral line, differ but not considerably from one model to other. 
However, these two models give differences on the profile of strong lines like A 4045 A 
(see Fig. 2). The enhancement factor of the 7h van der Waals constant has to be 2.5 for 
Holweger-Miiller's model and only 1.3 when one uses the Gustafsson's model, to fit perfectly 
the profiles by Kurucz et al. (1984). The value of 2.5 is classical and universally adopted for 
most of the lines, but it can vary from one line to another. This problem has been treated 
by Anstee, O'Mara & Ross (1997) who reproduce well enhancement factors ranging from 
1.4 to 3.3. Once that for late-type stars analysis we used the Bell et al.'s grid of models, we 
decided to keep the value of 1.3 for all radiative transition lines of our Fe I and Fe II models, 
since we do not try to reproduce perfectly all the line profiles. This problem was also found 
for the Ca I triplet fines (e.g. Cayrel et al. 1996), where Gustafsson's model looks less good 
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when using a theoretical enhancement factor of 2.44 for the Ca I hne A 6162 A, compared 
to the resuh obtained with the Holweger-Miiller's modeL This of course has an important 
consequence on detailed analysis of very metal-poor stars, because these strong lines would 
have to be used once that they are the only ones measurable on observed spectra. 

For solar photosphere, there were no important NLTE effects (overionization around 
0.02 dex) found by the position of the curves of growth of Fe I and Fe II. Once the 
precision of curve of growth's fit is around 0.04 dex, we have adopted 0.0 dex NLTE effect 
for Gustafsson's solar model. Consequently, one can say that the ionization equilibrium 
is perfectly reproduced by the LTE surface gravity logg© = 4.44. This absence of 
overionisation is due to important UV line blocking in solar dwarf stars and is independent 
of the solar model used. Fig. 3 shows the solar theoretical curves of growth of Fe I and 
Fe II considering only computed lines ranging between AA 2200 - 10000 A (see section 2.3). 
As we can see, curves of growth have a thickness which is produced by small NLTE effects 
giving an increase or a decrease of W, depending on population of levels from which lines 
are formed - e.g. depending on the excitation potential (see also Holweger 1996). These 
effects can be seen by the analisys of each line profile, whose detailed study is not the 
purpose of this paper, as mentioned in section 2.1. Our interest is centered in the problem 
of overionization, which have direct consequences on the ionization equilibrium (not in 
NLTE effects in line profiles) and therefore on curve of growth analysis. 

In 1949, Carter was the first to point out that the solar curve of growth in its damping 
part is divided into two main branches. Different conclusions were proposed to interpret it: 
the odd-even effect, the dependence of the damping constant on the multiplet of the lines 
(Carter 1949, Cayrel de Strobel 1966, Pagel 1965), but no comprehensible solutions were 
demonstrated (Foy 1972). Rutten & Zwaan (1983) had suspected NLTE effects. Usually 
neglected, this effect had no disastrous consequences on curve of growth analysis, because 
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the lines used to determine stellar abundances were not strong enough to lie in the damping 
part where this effect exists. Fig. 3 shows clearly double branches in the damping part for 
Fe I and Fe II. We remarked that the lower branch of damping curve of, for example, Fe I 
refers to lines originated from the ground levels (5D). The upper part is populated by lines 
originated from the (5F), which can be considered as resonance levels, because there are no 
permitted transitions with the ground level (5D). It seems that these 5D levels are pumped 
toward upper levels with more efficiency than the next levels (5F). UV lines clearly pump 
these levels as shown on Fig. 4, where bi = nf'""'"^/nj""'"^ is the classical coefficient used to 
study NLTE effect. In these two levels, bi have respectively values close to 0.13 and 0.24. 
Damping lines vary as the square of abundance, which means that for two lines having the 
same abscissae (log X) the ratio of their Wnlte is varying as the square of bj/bj, i.e. 0.14 
as can be checked on Fig. 3. This amplitude of the double branching of the Carter's effect 
is of the same order of magnitude as found on empirical solar iron curve of growth. We 
drawn attention to spectroscopists using lines originated from these Fe I resonance levels to 
be careful before deriving stellar iron abundances for very metal-poor stars, for which these 
lines are the only one easily measurable. Fe II lines present the same effects of splitting, 
but NLTE effects (see Fig. 4) are globally less pronounced, because they are formed in 
deeper parts of the photosphere. However, Anstee, O'Mara & Ross (1997) proposed a new 
approach of the theory of damping constant which could help to understand Carter's effect. 
Probably more detailed and precise computations are needed to determine the contribution 
of both results in Carter's effect. 



4. Comparison with other works 

Before to compute numerous new NLTE surface gravities logg and [Fe/H] for 
metal-poor stars, we decided to compare our technique predictions for two stars previously 
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studied by other authors. One is Arcturus (Takeda 1991) and the other is Vega (Gigas 
1986). In the case of the cool giant metal-poor star Arcturus we found negligible NETL 
effects of the same order as computed by Takeda: 0.03 and -0.02 respectively for Fe I 
and Fe II. For the hot star Vega, significant NLTE effects were derived. Because NLTE 
corrections on abundances of Fe I and Fe II were not the same, ionization equilibrium was 
not satisfied by the input of LTE surface gravity. We had to iterate until reach the following 
values for Vega : loggvega = 4.29, [Fe/H]vega = —0.33 (A0 = 7.46). This iron abundance is 
comparable to that of Gigas (1986), who proposed -0.55 dex (using a solar reference value 
of 7.67). It should be noted that most of Fe I lines are formed around logr ~ —1. — 2., on 
the contrary Fe II lines are generally formed in deeper parts of atmospheres. 

5. Results on metal-poor subgiant to subdwarf stars and consequences 

We selected in Thevenin's catalogue a set of 136 subgiant to subdwarf stars, with 
abundances ranging between -4.0 and 0.0 dex. We applied our strategy described in section 
2.3 to each of them. Results are presented in table 1. As one check of our final results, 
we compare computed Wnlte for Fe I and Fe II lines for the metal poor star HD 140283 
with those measured by Ryan, Norris & Bessel (1991). The correlation is shown on Fig. 
7. One can note that Fe I and Fe II lines are mixed, meaning that the estimated value 
logg = 3.74 is correct - same conclusion for [Fe/H] = -2.21. Its surface gravity was changed 
from giant to subgiant as derived by Hipparcos parallax (logg = 3.79, Nissen, H0g & 
Schuster (1997)). We show on Fig. 5 and Fig. 6 the amplitude of the overionization in the 
atmosphere of HD 140283 analyzed with the LTE atmospheric parameters. 

Nissen et al. (1997) have shown that exists a discrepancy between spectroscopic logg 
taken in the litterature and those deduced from Hipparcos parallaxes. On Fig. 8 are plotted 
the loggNLTE values derived by us versus logg Hipparcos (Nissen et al. 1997 and Clementini 



-14- 



et al. 1998) for stars in common, showing that our resuhs are remarkably close to those 
from Hipparcos. Among the four points having the worst correlation (Nissen et aL's values), 
two stars are suspected to be double. The adopted error bars are those of Nissen et al. 
(1997), Clementini et al. (1998) and Thevenin (1998). 

Overionization by UV lines seems to play an important role in stellar atmospheres of 
poor Late- Type stars. As can be seen on Fig. 6, bi Fe I coefficients are for most of the levels 
far below 1.0, increasing from resonance levels to upper levels until ~ 1.2 - upper levels are 
concerned by infra-red transitions. In consequence, the source function S^, ~ bj/b; x B^, 
satisfies the relation Sj, > B^, for most of UV transitions (where B^, is the Planck function). 
Also for strong resonance lines the mean intensity satisfies the relation J^, > Bj, and drain 
lower levels towards upper levels, which are more easily ionized. Consequently, upper 
levels are overpopulated, Jj, decreases far below the Planck function and produces infrared 
recombinations. These important mechanisms are well explained in Bruls, Rutten & 
Shchukina (1992) and we refer to this paper for more details. 

The results for our sample of stars are : 

- If Tcfj increases, UV radiation field and the pumping of the resonance levels increase, 
but the infrared radiation field decreases too and the overionization do not increase as it 
would be expected. 

- If Tcff decreases, infrared radiation increases and are more efficient to overionize the 
upper levels. The effect of Teg variations on the overionization for F to G stars is found to 
be not very pronounced. 

- If abundance decreases, UV line blocking decreases rapidly and the overionization 
becomes very important as shown on the relation between metallicity correction factors 
^[Fc/H] and metallicities estimated by LTE approach (Fig. 9). This reveals that the 
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overionisation increases rapidly when the abundance decreases from ^ —0.3 to -1.5 dex 
to reach a maximum for metal-poor stars having [Fe/H] ^ -3.0. There are no significant 
overionization for solar type stars. 

The most important parameter to produce overionization, the main source of NLTE 
effects on curves of growth technique, is therefore the variation of metal abundance in 
atmospheres of Late- Type stars. 

The precision of [Fel/H] or [Fell/H] determination by using theoretical curves of 
growth is estimated to be 0.04 dex, as mencioned in section 3, meaning that corrections on 
logg reach a precision of 0.08 dex. 

We computed Fe I and Fe II curves of growth for hot supergiant star parameters 
corresponding to stars analysed in the Magellanic Clouds. These stars have moderate 
underabundances (~ —0.6 to —0.2) dex. We found a [Fel/Fell] balance between 0.05 
and 0.01 dex for SMC stars (see stellar parameters in Thevenin 1998). Holweger (1990) 
mentioned that iron NLTE effects computed using the LTE atmospheric parameters of 
the poor giant star CD — 38°245 lead to an error in the ionization equilibrium balance of 
+0.3 dex, value very similar to those derived by us for dwarfs stars. These results are not 
surprising as mentioned by Feltzing & Gustafsson (1998), once the UV opacity increases 
with the increase of the Balmer lines when the surface gravity decreases. This only means 
that the idea that NLTE deviations must necessarily increase with the decrease of surface 
gravities (hence increase in collisions) is not entirely correct. From the point of view of 
Boltzmann's law this is true, but if the UV flux is blocked by increasing of Balmer lines, 
deviation from the point of view of Saha's law do not increase with logg decreasing. 
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6. Conclusion 

We presented a study of departures from LTE for Fe I and Fe II, mostly for metal-poor 
stars. The mechanism of these departure is clearly identified as overionisation, responsible 
for important corrections on the values of surface gravities. A good correlation between 
our derived surface gravities and those deduced from Hipparcos parallaxes is a proof of 
the validity of our results in table 1, having important consequences on distances in the 
Galaxy. These corrections on the surface gravity could have also important incidences 
on the abundances of elements like Be (Gilmore, Edvardsson & Nissen 1991). NLTE 
abundance corrections are found to be less than 0.35 dex, but not negligible. Therefore, 
stellar abundance ratios could have to be revisited after having estimated possible NLTE 
effect on other elements like Ca, Mg, Al, O and Be. We are preparing papers on these 
subjects. We recommend to stellar spectroscopists working on metal-poor stars to use 
NLTE computations or the table 1 before to publish their LTE results, or to use surface 
gravities derived by Hipparcos parallaxes combined with a LTE analysis of Fe II lines which 
do not suffer important NLTE effects. 
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Fig. 1. — Grotrian diagrams for Fe I and Fe II atoms. Each configuration term is denoted 
by an alphabetic notation, representing the degeneracy of the state. Two Fe II levels do not 
have configuration terms identified, then they were labeled 62065 and 65363 (which refers to 
line transitions AA 6206.5 and 6536.3 A) as Hirata & Horaguchi's table. 

Fig. 2. — Comparison between two calculated Fe I NLTE line profiles and observed spectra for 
the sun. The theoretical profiles were calculated using the atmospheric models of Holweger- 
Miiller (1974) and Gustafsson (1981), for different enhancement factors j'h of the van der 
Waals constant. Observed spectra are from Kurucz et al. (1984). 

Fig. 3. — Theoretical Fe I and Fe II curves of growth for the sun. logX = loggf + logF + 
log(NFe/NH), where gf is the oscillator strength; F is a parameter defined for each line 
as a function of the atmospheric model (Cayrel & Jugaku 1963); Npe/Nn is the relative 
abundance. 

Fig. 4. — Departure coefficients of Fe I and Fe II versus optical depth in the solar case, 
bi = n|^^"'"^/nj""'"^, where nf'""'"^ and nY^^ are NLTE and LTE populations respectively for 
each level. 

Fig. 5. — Theoretical Fe I and Fe II curves of growth for the star HD 140283, computed 
with LTE atmospheric parameters: ^cff=0.90, logg = 3.20 and [Fe/H]=-2.5. For Fe I is 
overplotted the classical LTE curve of growth (full line), showing the need of correction by 
NTLE effects. 

Fig. 6. — Departure coefficients of Fe I versus optical depth for HD 140283 (^eflf=0.90, 
logg = 3.20 and [Fe/H]=-2.5). h-, = nf^^'^/n^™, where nf^TE ^nd n^™ ^^^ j^lTE and LTE 
populations respectively for each level. 

Fig. 7. — Computed equivalent widths for Fe I and Fe II lines versus observed ones by Ryan 
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et al. (1991) for Fe I and Fe II for HD 140283. Atmospheric parameters for calculation of 
equivalent widths are corrected of NLTE effects. 

Fig. 8. — Comparison between our derived log g and those using Hipparcos distances from 
Nissen et al. (1997) and Clementini et al. (1998). loggNLTE error bars are estimated from 
classical detailed analysis for dwarf stars ~ 0.30 dex (Thevenin 1998). 

Fig. 9. — Final estimated NLTE abundance corrections Ajpe/H] versus LTE atmospheric 
parameters from Thevenin (1998). 
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